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Abstract 

We prove the unique existence of solutions of the 3D incompressible Navier-Stokes equations 
in an exterior domain with small non-decaying boundary data, for t £ M. or t G {0, oo). In the 
latter case it is coupled with small initial data in weak . As a corollary, the unique existence 
of time-periodic solutions is shown for the small periodic boundary data. We next show that 
the spatial asymptotics of the periodic solution is given by the same Landau solution at all 
times. Lastly we show that if the boundary datum is time-periodic and the initial datum is 
asymptotically discretely self-similar, then the solution is asymptotically the sum of a time- 
periodic vector field and a forward discretely self-similar vector field as time goes to infinity. It 
in particular shows the stability of periodic solutions in a local sense. 

Keywords: Navier-Stokes equations, exterior domain, spatial asymptotics, time asymptotics, 
time-periodic. Landau solution, discretely self-similar, stability. 
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1 Introduction 

Let C K'^ be an exterior domain with smooth boundary 957 and unit outernormal N , and / = R 
or J = (0,00) be the time interval. In / x 51 we consider the nonstationary Navier-Stokes equations 

5tu - Am + (7i • V)ii + Vp /, divu = 0, (1.1) 

u\dn = u.^,, lim it = 0, (1.2) 

I >oo 

where uit, x) : / x — > M'' is the unknown velocity field, p{t, x) : / x — > R the unknown pressure, 
the given boundary data and f = fo + y ■ F the given force with F ~ [Fij ) being a 2-tensor and 
(V • F)j = diFij. We use summation convention for repeated index and u • V = In the case 

/ = (0, 00), we add the initial condition 

M|t=o Wo, divuo = 0. (1.3) 

By an exterior domain we mean a connected open set with bounded complement. Without loss of 
generality, we may assume the complement of H contains the origin and is a subset of B]^^ := {x g 
R3 : |a;| < for some i?i > 0. 

In this paper, we first consider the solvability of the problem (jl.l[) - (|1.3|) with nondecaying bound- 
ary data in time. Moreover, we also study the asymptotic properties of solutions with time-periodic 
boundary data. 

To explain the background, we start with the review of the stationary problem. In 1965 Finn [9] 
showed the existence of a small stationary solution satisfying |u(a;)| < C|a;|~^ in £7 for small data. 
Nazarov and Pileckas [57] proved that the solution is asymptotically self-similar at spatial infinity, 
i.e., the solution converges to a (-1) -homogeneous vector field faster than C|a;|^^. Recently Korolev 
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and Sverak [17j showed the asymptotic profile is given by a Landau solution. In particular the decay 
rate \x\~^ is optimal in general. 

Landau solutions are a family of vector fields U'^ and functions in , with parameter 6 € M'^ , 
that solve 

- Au + (m • V)u + Vp = 6<5o, div u = 0, (1.4) 

where 5q is the delta function at the origin. They are axisymmetric: in spherical coordinates p, 9, ip 
with b in the direction of the north pole, 



^ / - 1 \ 2 sin 03 , ^ 

^ =-[t-a ^-l^P + Oee-^^ ^e^, 1.5 



where A = A{\b\) e (l,oo] is determined by |6| = IQt:{A + \A^\og^^ + -^(jtzT)) and strictly 
decreasing in \b\ > 0. It is the unique solution of (|1.4p in the class of (— l)-homogencous vector fields 
in R3. See Landau HHHS], Tian and Xin [55], and Sverak [M] . 

On the other hand, there are also a lot of works on the time-periodic solutions in exterior domains. 
After the earlier works by Salvi [3D] and Maremonti-Padula [23], Yamazaki [37] showed the unique 
existence of time-periodic solutions in the Lorentz space L^^°° for zero boundary data and small 
forces F G 7^3/2,00 -^^rj^j^ — Galdi and Sohr [12] further showed the existence of time-periodic 
solutions satisfying the pointwise estimate |M(t, a;)| < C|a::|~^ where C is independent of time. 

When we are interested in time-periodic solutions or, more generally, the exterior problem for 
nondecaying boundary data and forces, the function spaces should allow nondecaying functions in 
time. On the other hand, in view of the optimal decay rate for the stationary solutions, it is 

natural to choose the spaces L^'°° (weak-i^ space) or Xi in spatial variables, where Xk (for A: > 0) 
is the space of functions defined by the norm 

II^IU, :=sup(l + |x|)'=Kx)|. (1.6) 

Let BCw{J] X) be the class of bounded and weak-star continuous X-valued functions defined on a 
time interval J. We omit the subscript if it is strongly continuous. 

Our first result concerns the unique existence of very weak solutions, to be defined in ^32. 21 with 
small nondecaying boundary data in the time interval M or (0, 00). As a consequence, we obtain an 
existence theorem of time-periodic solutions. 

Theorem 1.1 (Existence) Let fl be a smooth exterior domain in with d^l C {x : \x\ < Ri}. 

Let / fee M or (0, 00). Let S be a positive constant. There are constants Eq > Eq > and C > such 
that the following holds: Let uq, u^, /q and F be given data with the convention uq ~ if I — M,. 
(i) Assume that 

e := ||u*||n/i,o=(/;C2(an)) + ll/o||L~(7;X3+i) + \\F\\l^(i-l3/'^.'^{q.)) + \\uo\\Li--^(n) < £o- (1-7) 

Then there is a unique very weak solution u E BCw{I', L'^'°° (fl)) of (jl.ip - (|1.3p satisfying 

ll'«||L~(/;L3.°°(n)) < Ce. (1.8) 

(a) Let / = M. Assume that all u^,, /o and F are periodic in time with period T > 0, then the 
solution in (i) is also periodic in time with same period. 
(Hi) Let / = M. Assume that 

e W2.==(R;C2(ao)) + II/o||h'1.==(R;X3+s) + ll-P^II Wi.~(R;X2) < ^0, (1-9) 

then the solution in (i) satisfies 

\u{t,x)\<Ce\x\-^, (|x| > t e K). (1.10) 
Comments for Theorem lj.il 



1. The initial condition (|1.3p is understood in the weak-star sense. Thus we do not need a 
compatibility condition between u^^ and uq. When u^^ = 0, we usually also require u belongs 
Ll'°°{D) which is the subspace of L^'^{Vt;m.^) with divu = and u ■ N\an ^ 0, N being 
the outer-normal of fl, see ^2.11 This is not suitable for the inhomogeneous boundary value 
problem. 
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2. In the case J = M and is time-independent or timc-pcriodic, wc recover the rcsuhs of Finn 
[9], Galdi and Sohr [12] on the existence of solutions behaving hke |a;|~^ as |a;| oo. 

3. As noted above, there are a lot of literature on the existence of periodic exterior flows for 
zero boundary data [301 [24l [37l [12] . Our result allows nonzero boundary data. There are also 
results on existence of solutions for nonzero boundary data e.g. Amann [2], Farwig, Kozono 
and Sohr [^. However, most of them require time decay of the boundary data, which is not 
suitable for the periodic solutions. Since the maximal regularity estimates in [5] [5] arc not 
available in our solution spaces, we use the duality argument by Yamazaki [37] to construct 
the solution. In fact, we first decompose u = E + v where E is an extension of the boundary 
data u*, and then wc construct the unique mild solution v of the difference equation (|3.8p 
following the method by [37] • Unlike [37], we do not require time continuity of the force (in 
part (i)) and our solution is only weak-star continuous in t. We then show the equivalence of 
this solution and a very weak solution in ^13. 21 Note the datum u^, is in the class, not usual 
C^'", < a < 1. The existence of the extension E is shown in Lemma [131 

Our second result concerns the spatial asymptotics of time-periodic solutions. To describe it, we 
recall the momentum flux density tensor for a solution {u,p) of (|l.ip . 

Tij{u,p,F) =^ pSij + UiUj - diUj - djU^ ~ Fij. (l-H) 

Equation (|l.ip can be written as 

dtuj + d^j = /oj , div u = 0. (1.12) 

We show that the asymptotic profile of a time-periodic solution is given by a Landau solution 
determined by the tensor and chosen independent of time. 

Theorem 1.2 (Spatial asymptotics of time-periodic solutions) For any T > 0, R > and 

a G (1,2), there are constants ei > and C > such that the following holds. Suppose {u,p) is a 
time-periodic solution of (|l.ip with period T for i? < |x| < cxj, and satisfies 



£:=sup sup {|x|2+"|/o(i,a;)|-f |a;^+"|i^(^,a;)|} + h(^,•)||L3.-(|,|>fl)^- 
*.«lH>i. ^^_^3^ 

+ sup {\u{t,x)\ + \Vu{t,x)\ + \p{t,x)\]\ <£l. 
FI.<\x\<R.+ l ) 

Then \u(t,x)\ < C£|a;|~^ for \x\ > R and all i e R. Moreover, let Tij be defined by (|l.lip . let 

6, = lim i / / T.judS^dt, (n, = ^, 1-1,2,3), (1.14) 



Ja J\x\=p F. 
and let be the Landau solution corresponding to b, whose choice is independent oft. Then 

\u{t,x)~U''{x)\<Ce\x\-"', {\x\> R; t gR). (1.15) 
Comments for Theorem 

1. In Theorem 11.21 we do not specify a boundary condition for u. All we need is a solution for 

> R. Thus our result is applicable to the periodic solutions constructed in Theorem 11.11 
(ii), (iii). 

2. If u is independent of time. Theorem 11.21 recovers the result of Korolev-Sverak [T7| for small 
exterior stationary Navier-Stokes flows. We observe that the decay rate |u(i,a;)| < C|a;|"^ for 
small periodic solutions is optimal as well as the stationary case. 

3. In order to see that the limit in (|1.14p exists, denote the integral as Ij{p) and then we have 
from ([02]) that 

IM - Ijipi) r [ fo, {R<Pi< P2). (1.16) 

^ Jo Jpi<\x\<p2 

We would like to emphasize that the choice of b in (|1.15p is independent of t. There is some 
sort of cancellation effect behind it. 
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Finally we consider the large time asymptotics of solutions of the problem p.ip - (jl.3p when the 
boundary data and the force are time-periodic. Borchers and Miyakawa [3] showed the stability 
of stationary exterior flows under small initial perturbation in L^^°° vanishing at boundary. For 
earlier stability results, see e.g. [Tallin]. Our third theorem extends these stability results to small 
time-periodic solutions. 

Theorem 1.3 (Time asymptotics) For any r>0, 5>0, 77>0 and 3 < gi < there is 

£2 > such that the following holds. Let and f be time-periodic data satisfying (|1.9p . Assume 
the initial data uq is asymptotically self-similar in the sense that there is a (- 1) -homogeneous vector 
field uq e Xi such that div uq ~ and 

lluollxi < £2, Vuo e X2 and e -.^ \\uq - uo\\ ^ 3 <e2- (1-17) 

Then the solution u in Theorem \l.l\ (i) for t > can be decomposed as 

u = Q + w + r. (1.18) 

Here Q is the periodic solution for the data and f in Theorem \l.l\ (Hi). The term w is the forward 
self-similar solution of the perturbed Navier-Stokes system in (0,oo) x R'^, 

9tw - Aw -f V(w «) w H- «) w + w «) C/^) + Vpo = 0, divu; = 0, (1.19) 

with initial data wq ~ uq ~ , where U''{x) is the Landau solution corresponding to Q given in 
Theorem M.'A with a = 1 + 5 , and w satisfies 

\w{t, x)\ < Ce2{\x\ + ^/^)-l+''|:E^^ {t e (0, 00); x e R^). 

The term r satisfies the following decay estimate: 

IkWllL.-(o) <Cei-^(*-i)^5 vt>0, Vge[3,gi]. 

Comments for Theorem \1.3\: 

1. The last assumption in p.l7p means that the profile of uq at spatial infinity is given by uq. 
Since ||(5(t, a;)||L9.«= < C, \\w{t,x)\\L'i.^ < Ct^"^!-?^^) for 3 < q < qi and r decays faster than 
the other terms as t — > 00, our theorem shows that u is asymptotically equal to the sum of Q 
and w. In particular, this implies the stability of the periodic solution. 

2. In the case f2 = with zero force and no boundary data, Planchon P2] showed that if the 
initial data is asymptotically self-similar, the solution is also asymptotically self-similar. Here 
the asymptotic profile of the solution at large time is given by the self-similar solution of the 
non-perturbed equations (jl.ip with zero force. See also [?]■ In our case, the asymptotic profile 
w satisfies an equation modified by the Landau solution but it is still self-similar. Unique 
existence of the self-similar solutions for (|1.19p is considered in Proposition l5.ll 

3. If Uq is not self-similar, we still have similar decomposition like (|1.19p . In particular, if uq is 
asymptotically discretely self-similar, w becomes forward discretely self-similar. The notion of 
the discretely self-similar solution is discussed in the introduction of section [5j In section [51 
we will show a more general decomposition. See Theorem 16. II 

Our previous three theorems can be similarly posed in the entire with a singular force sup- 
ported at the origin 

^tw- Aw + (ii- V)w + Vp = 6(t)5o, divw = 0, (a; e R^), (1.20) 

which is studied by Cannone and Karch [J . The problem will be pursued elsewhere. 

This paper is organized as follows. In section [2J we recall some preliminary results related to the 
exterior problem. In section [3] we prove Theorem 11.11 Section 2] is devoted to the proof of Theorem 
11.21 In section [5] we consider the perturbed Navier-Stokes system (jl.l9p which plays an important 
role in the next section. Finally we consider the time asymptotics of the exterior fiows and prove 
Theorem ll.3l in section [51 

Notation, a < b means a < Cb for some constant C . a ^ b means a ^ b ^ a. (a) = -^/jap + 1. 
For 1 < p < 00, its conjugate exponent p' is defined by 1/p + 1/p' — 1. We denote M+ = (0, 00). 
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2 Preliminaries 

In this preparation section we recall Helmholtz decomposition and Stokes semigroup in ^2.11 define 
very weak solutions in ^2.21 prove extension lemmas in show that /o can be absorbed into VF 
in H'lAl and prove decay estimates for Stokes system in ^2.51 

Recall rt denotes an exterior domain in M.^ with smooth boundary dil and unit outer-normal N. 
We also assume ^ Ct and dfl C Bji-^ . 

2.1 Helmholtz decomposition and Stokes operator 

The Helmholtz decomposition 

L'?(17;R") = L«(f7)eG''(f}), (l<q<oo), (2.1) 

for a C-'^-exterior domain f2 C M", n > 2, is well understood. It is first proved for n = 3 by Miyakawa 
[26] and for all n > 2 by Simader-Sohr [31]. Let P = Pg be the associated Helmholtz projector from 
L'^ onto L^. Then P can be extended by interpolation to a bounded projector on each Lorentz space 
L'^-^{n), 1 < q < oo, 1 < r < oo, with the Helmholtz decomposition 

L«''^(17;M") = L9''-(f])©G«'''(f}), (2.2) 

where 



(2.3) 



Ll'''{n) ^{ue L'i ''{n- R") : divu = 0, u • N\an = 0}, 
G^^^in) = {Vp e L«''^(0;M") : p € Lf^.m. 

Furthermore, if 1 < r < oo, then {LY)* = l'-/'"'"^ and (G'?''')* = G(«'''''). See [1 Th. 5.2]. 
The Stokes operator Aq on L^{^) with the dense domain 

v{Aq) = Li{n) n Wo''^{n) n w^'^n) (2.4) 

is defined by AgU — —PgAu for u e D{Aq). It extends to a closed linear operator on LY{^l) with 
domain 

V{Aq^r) = {ue Ll'^m : V^u e L«''^(n), J = l,2,u\on = 0}. (2.5) 

One also has 

'D{Ai/;) = Lrmnw^^^'^^'-\n). (2.6) 

The semigroup {e~*^}t>o also extends to i^''"(il). 
Wc now recall some estimates in Lorentz spaces. 

Lemma 2.1 Let fl gM.^ be an exterior domain with smooth boundary. One has 

WA^^^uWLP.r <\\Vu\\Lp.r, (Kp < c5o,l < r < oo), (2.7) 

I|Vu||lp- < (l<p<3), (2.8) 

||V«|Up,i < (l<p<3). (2.9) 

The above is [2Z1 Theorem 2.1]. 

For 1 < p < 3, we define p* by 1/p* = 1/p - 1/3. 

Lemma 2.2 Let fl CZM.^ be an exterior domain with smooth boundary. 

(i) For l<p<q<oo,a^ |(i - i) and 1 <r < oo, for any (j) € P^''"(fi) 

l|e"*^0llL.-- < i""!l<^||Lp-, t>0, (2.10) 

< t-"-^||(/)|Up.^, OO. (2.11) 
(a) For 1 < p < 3, for some c — c($l,p) > 0, for any cj) G Lg^^(il), 

/ \\A^/'e-'Uh,^.^dt<cUhp... (2.12) 
Jo 
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(Hi) For 1 < p < 3 and < a < 1, for some c ~ c(fl,p,a) > 0, for any (f) G LP'^{fl), 

-^\\A^^^e-^*-'^^(l)\\LP',ids<ct-^U\\L,,i, t>0. (2.13) 



t 

s 



If^<P<%,'we may replace A^l"^ in ((2?T2|) and ((2?T3)) hy V due to ^3, 



Proof. The decay estimates (^1^ amd ([^IT]) are [371 (2.4), (2.5)]. (Although [37] only states 
the case r = 1, it works for other r.) The integral estimate (|2.12p follows from [37l (2.11)]. 

To show (PT3)) . divide the integral to s < t/2 and s > t/2. Applying (pHI) for s < t/2 and 
([212)1 for s > t/2, we get ([21^ . □ 

As a corollary, we have the following lemma. Here BCw denotes the class of bounded and 
weak-star continuous functions. 

Lemma 2.3 Let fl CZMp be an exterior domain with smooth boundary. Let / = M or L = M_(_. For 
fit) e L°°(/;L9'°°(r2)), 3/2 < 9 < 3, define 

{Gifm = f e-^'-^^^PWf{s)ds, [t e L) (2.14) 

Jinf / 

in the sense that, for p = (q*)' , that is, 1/p — 4/3 — 1/q, 

mfm, 0) = r (/(s), Ve-(*-^)^0)ds, V0 € LP/, yt e /. (2.15) 

Jinfl 

Then Gif e BC.u,{i;Li-°°{VL)), and for some c = c{n,q) > 

ll^//llLoo(/;Lr-) <c||/I|l~(/;L..==>)- (2.16) 

Note the integral Gif is in weak sense and may not converge absolutely. 
Proof. For any e L^'^ and t G I, we have 



t-inf / 



= l / (/(t-r),Ve--^0)dr| 

< / \\f{t-T)\\L.-^\\We-^^^^,>.^dT. 

Jo 



(2.17) 



Note that 3/2 < q < 3 is equivalent to 1 < p < 3/2 and that q' ~ p*. Then Lemma [221 (ii) implies 
that 

1(57/(0,0)1 <ci|/||L-(/;L.-)li'/'liLp.i, (2.18) 

which shows Gif G L°°{L; L% '°°(ri)) and (|2.16p . To show weak continuity, a computation similar to 
the above shows Vt < t' E I and e := t' — t 

\{Gif{t')-Gif{t),cj,)\ 

<J^ \{f{t'-r),ye-^^cj^)\dr + J^ |(/(t - r), Ve--^(e-^ - l)0)Mr ^^.19) 



<ll/l|™-^ ||Ve--^0||^,,,,dr + ||Ve'-^(e-^'^-l)0|L,,.,dr 

By Lemma [2IH and by the strong continuity at e = of the Stokes semigroup in L^'^, the right hand 
side converges to as e — )■ 0+, and hence \{Gif{t') — Gif{t), ^ with cither t or t' fixed. □ 

2.2 Very weak solutions 

In this subsection we define very weak solutions in an exterior domain Q with unit outernormal N. 
Our definition is a variant of that in Amann [2] and Farwig, Kozono and Sohr 0. We take the 
constant k to be 1 for the Navier-Stokes system, and k = for the Stokes system. 
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For finite or infinite time interval J, denote tlie space of test functions 



V,:=\v.^Cinr, R-) : ^^^.^ ^ \^ '^'^ A , (2.20) 

I ^ ^ divw = 0, suppwdOj, w\on=0<' ^ ' 

wliere flj ^ J x Q contains tlie boundary. Tlius Vw may be nonzero on boundary, and w(to, •) may 
be nonzero if to = inf J is finite. 

Let / = M or / = R+. A vector field u G Lf^^{I x 17; M^) ig called a very weak solution of the 
Navier-Stokes system (when k = 1) or of the Stokes system (when k = 0) with initial datum uq (with 
the convention uq = if / = K), boundary datum m*, force / and mass source k (with sufficient 
regularity) if 



{ — (it , wt + Aw)j^ + {u^, , N ■ Vw)gj^ — {ku u , V?«)fj — {ku , w)^^} dt 
= (wo , w'(0))o + / (/ , w),^dt, e P/, 



(2.21) 



and 



divu(t) = fc(t), TV • u(t)|ao = iV • it*(t) for a.e. < e /. (2.22) 
An elementary calculation shows that for w G 2?j 

TV • Vw = curlw X TV on (2.23) 

Thus (|2.2ip contains a condition only on the tangential component Nxu^ of on dQ, and we have to 
assume the additional condition in (|2.22p for the normal component N-u\gi-i = N-u^,. Note that, when 
O is a bounded domain, one needs to assume the compatibility condition Jg^ (t) ■ NdS = k(t)dx 
for a.e. t. When D, is an exterior domain, however, this is unnecessary. For the rest of this paper we 
take fc = 0. 



2.3 Extension lemmas 

We will use the following extension lemma to extend a given boundary data m* with zero flux on 
every connected component of dfl to a divergence-free vector field of compact support in f2. The 
special case we need is u* € C^{dil) and J7 is a smooth exterior domain in R'^. 

Lemma 2.4 Let I G N and < a < 1. Assume is a domain in M", n > 2, with compact 
boundary 9il of class (7'+^'" and unit outernormal N . Suppose dil has M connected components 
Tk, fc = 1, . . . , M . For any 5 > Q, there is a linear map £ which assigns for each G C''"(9f2, M") 
with Jp^ Uif ■ N — Q for all k, a vector field u — £(7i*) G C''"(j7) so that 

divu = 0, u|ao=w*, suppuC^s, \\u\\ci.^{n) < C\\u^\\ci.<^(dn)- (2.24) 

Above = {x G ri : dist(x, Oi7) < 5} and C ~ C(il,l,a,S). This linear map restricted to smooth 
is the same for all {I, a) so that dil G C'^^'". 

The following is a more general result which we will not use. 

Lemma 2.5 Let I G N and < a < 1. Assume n is a domain in M", n > 2, with compact boundary 
d^l of class (7'+^'" and unit outernormal N . For any open convex set uj C M" containing dil, there 
is a linear map £ which assigns for each G C''"(9ri, R") with Jg^u^, ■ N ~ 0, a vector field 
u = £{u^) G C''"(17) so that 

divu = 0, suppu C ClDuj, ||w|lc''°(o) < C'll"*llc''°(ao)- (2.25) 

This linear map restricted to smooth is the same for all {I, a) so that d^l G C'"*"^'". 

Remarks on Lemmas \2.4\ and \2.5\ 
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1. If G C^'"(5rJ), < a < 1, Lemmas 12.41 and 12.51 fLcmnia [^751 in the case ft is bounded) 
is proved in Kapitanskii-Pilctskas |15| . In addition to induction in dimension, they write 
Ui = di{wij + 6ij(f>) with Wij anti-symmetric and estimate Newtonian potentials in Holder 
spaces. This is why they do not allow a = or a = 1, which are allowed in Lemmas 12.41 and 



2. In case 51 is an exterior domain, the extension constructed in |15| does not have compact 
support unless one further assumes /p^ • = for all k. See also Kozono-Yanagisawa [H] 
for extensions in more general domains. 

3. A related problem is the construction of a vector field v for a given / : 51 — > M satisfying 

/ = so that 

divv = f, {xen), v\an^O, ||Vi;|| < C||/|| (2.26) 

for suitable norms. The extension problem is reduced to (|2.26|) as follows: first extend to 
U G C'°'"(ri) which may not be divergence- free, then solve (|2.26p with / = divU, and finally 
define u ^ U — v. The problem (|2.26p is solved first for / G C^{fl) and then for general / by 
approximation. For 1 < q < oo, this approach is good for it* G W^~^^''''^{dfl), but not suitable 
for G l^'=-i/9'9(ar2), k>2, since C^in) is not dense in W''^'^ f] Wq''^. Similarly it is not 
suitable for G C^{dfl). 

4. Remark VIII. 4.1 on [TTl p25] says that there is no linear map which assigns a v satisfying 
(|2.26|1 for a given / so that \\v\\li < c\\f\\iv-i,q, which we fully agree. However, we do not 
understand why it continues to assert that, "by the same token", there is no linear map which 
assigns a v for / G 0(^1) so that ||w||ci < c||/||c- 

5. Theorem 4 of [l5] asserts the existence a solution v of (|2.26p satisfying (for < a < 1) 

Il«llc^+i-(j2) <C||/||c^,»(n). (2.27) 

However, its proof assumes dist(supp/, dQ) > and the constant C depends on this distance. 
See [m §6]. 

Proof of Lemma \2.4\ We will prove the case n = 3 for notational simplicity. The general case 
(n > 2) is proved in the same way. Denote x' — {xi, X2) and the open square K = {x' G : — 1 < 

Xi,X2 < 1}. 

We first consider the case that the support of is on K x {0} with jK x (0, 1) C 51. In this case 
we have /^u^ = 0. Extend u^{x') ^ for x' G M.'^\K. Choose (j) G C°°(M.'^) so that (f){^) = 0(|^|), 
(/)(C) = for 1^1 < 1/8 and |^| > 1/4, and /^^ (^(O^C = 1- Also choose x e C°°(R) so that x(s) = 1 
for s < 1/8 and x(s) = for s > 1/4. We define the extension by 



U{x)^{d3^l, 93$2, -ai$l-a2$2), (2.28) 

wher^l 



<^j{x',X3) = xixs) 



, (.? = 1,2), (2.29) 



X3 / u^x' + X30md^ + Pi^') 
and, with g{xi) = Jg^u^{xi, X2)dx2: supported in < 1 and j\ g{xi)dxi = 0, 

f\x') = -X'ix2) I ' ff(ii)d5;i, f\x') = / ' ul{x^,x2)dx2 - g{xi){l - x{x2))- (2.30) 



(See the proof of Lemma [2.61 ) Note f^ are supported in K, div^a / = m^, and u is supported in 
I X [0, 1/4]. One verifies ^(TI^ directly. 

We next consider the case that the support of is on a graph over K: X3 = h{x') for x' £ K 
with \h{x')\ < 1/4, and 51 n {^K x [-1, 1]) lies in X3 > h{x'). We have 



TV = -(9i/i, (92/^, -1), e = a/|V/i|2 + 1, dA = edxidx2, (2.31) 



^The formula for is from Lemma 2.5. The formula for and «^ extends the n = 2 case of 1151 (7.3)]. 
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and = • NdA = Jj^u^, ■ (dih, d2h, ~l)dxidx2- Define new coordinates 

yi = xi, y2=X2, y3 = X3-h{x'), (2.32) 

and, for each vector field defined ioi x E Cl O {K x [0, 1]) define a new vector field for y G K x 
[-3/4,3/4]: 

U{y) = {u\x), u'^ix), u^{x)-u\x)dih{x')~u^{x)d2h{x')). (2.33) 

For a given defined on dft, we define U^, on K by the same formula. Then [/* satisfies Jj^U^ = 
and we can extend U from [/* by the previous case so that div?7 = and supp U C K x [—1/2, 1/2]. 
We finally define u from U by (|2.33p . One checks directly that div^; u = divj, J7 = 0. 

For the general case, we can find finitely many balls Bj, j = 1, . . . , J, with same radius p and 
concentric balls B* with double radius 2p, so that dil C ^jBj, LljB* C Bji^, and that dD. r\ B* is a 
graph in B* in direction /ij, and belongs to \ {x — Xj) ■ < p/8 where xj is the center of Bj. Choose 
a smooth partition of unity {r]j}j on dfl so that J2j Vj = 1 on dfl and suppT^j C Bj n dO,. 

For any given with /p • iV = for all k, we claim we can decompose 

u* = \] supp J C 917 n Bj, / • = 0, (2.34) 

j ' ' Jan 

with suitable estimates. When dft has only one component {M ~ 1), choose Xk, k ~ 1, . . . , J — 1, 
so that Xk S C^(9f2 Ci BkCi UjykBj) with J^j-^ x^- = 1. We define by induction: Let Ui ~ u*, 
and for j = l,...,J— 1, 

=Ujiij - Xj Ujijj, Uj+i=Uj-mj (2.35) 

and = C/j. When 9fi has more than one component, we can perform the above decomposition 
for each component. (The above decomposition follows the proof of [TUl Lemma IIL3.2].) 

Since Jg^ u^,j ■ N = 0, by the second case above with a suitable rescaling we can extend to 
divergence- free Uj supported in B* n Cls- We now define £{u^,) = J2j '^j- ^ 

Remark. If is prescribed in K x {0, 1} with u*(.t', 0) = (0, 0, vq{x')), u^{x' , 1) — (0, 0, vi{x')) 
and Jj^vo{x')dx' = J^^ vi{x')dx' , we can extend u^^ to K x [0, 1] by 

u{x) = ii;'ix3)f\x'), ij\x3)fix'), (1 - ^(x3))^;o(^') + ^(x3)«i(x')), (2.36) 

where tJj = 1 — Xi 

g{xi)dxi, f{x')^ I {vo ~ vi){xi,X2)dx2 - g{xi)'ilj{x2), (2.37) 

-OO J — OO 

and g{xi) = /jg(uo - ui)(.ti, s)ds. 

Proof of Lemma \2.5\ The case M = 1 follows from Lemma 12. 4[ with (5 > chosen so small 
that Vis C w. Suppose now M > 2. 

Claim: We can choose A/ — 1 line segments connecting Fj , so that intersect dft only at end 
points and at right angles, and each F^ intersects at least one Lk- This is chosen by induction. Let 
Ai = {1} and Am — {1, . . . , M}. Suppose Ak, a subset of Am has been chosen for k = 1, . . . , M — 1. 
Since d^l is compact, one can find a line segment Lk minimizing the distance 

dist(U,gA,F„ U,^A,r,). (2.38) 

Clearly Lk intersects Uj^A^^j only at one endpoint and Lk intersects Uj^A^Fj, say Fj^^^, at the 
other endpoint. Moreover, Lk without endpoints is inside fi. Now let Ak+i = Ak U {jk+i} and 
continue. 

Once all Lk have been chosen, we give a rank of F^ and Lk as follows: Tj is assigned rank one 
if it intersects only one Lk- This line segment is renumbered as Lj and also assigned rank 1. Let A 
be the set of Fj and Lk without those of rank one. Rank one boundaries and line segments in this 
reduced set are assigned rank 2. We continue this exercise until we exhaust all Fj and Lk- 
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Fix < £ ^ 1. For each k, let be the unit direction vector of Lk (unique up to a negative sign), 
let L|. be those points on the line extending with distance to Lk less than e, let rfc(x) = dist(a;, 
and let U^{x) = e^~"x(£~^''fe(2;))efe where x : K — > [0, 1] is a fixed smooth cut-ofF function with 
x{r) = 1 for r < 1 and x(^) = for r > 2. For e > sufficiently small (depending on dil), we have 
Jj^ U'^ ■ N = ±(Jjg„_i x{\x'\)dx' + o(l)) 7^ if Tj is one of the two boundaries intersecting Lk- We 

also have supp C w. 

Now, for a given it*, let u'^ ~ and define recursively for fc > 1: 



of rank k 



Denote the final one as . One verifies that Jp^ u!f • TV = for all k and + |cj| < 

C||u^^||ci,o. We now define 

A/ 

£:(u*) = f (uf ) + E CjC/^'(a;). (2-40) 

where £{u^^) is defined by Lemma \TM One verifies that ||f (M*)||c!,a(Q) < C||w^^j|pi,c« kjl — 

C||uf^||(7i,o, and that the support of £{u^) is inside w. □ 

2.4 Source terms 

In this subsection we show that any force f{x) in R" decaying like Ixl""""^ as infinity can be 
decomposed in the form / = /o + V • -F with supp /g being compact. 

Lemma 2.6 // f{x) is defined in R" with |/(a;)| < (x) °, a > n > 1, then for any R > we can 
rewrite 

n 

f{x) = fo{x)+J2d,F,{x) (2.41) 

where supp /o G Bi?(0) and \F,{x)\ < (x)-"+^||(.T)"/(a:)||L- • 

Remark. (i) If /(a;) is defined in an exterior domain C K" with ^ H, we may extend / 
by zero to entire R", and choose i? > so small that il n Bj^ = 0. Then the first term /o(a;) in the 
decomposition (|2.41l) can be ignored for x E ft. 

(ii) If / = go + V • G with |.go(a;)| ^ (x) a > n > 1, and |G(a;)| < (x) , we may assume go 
has compact support by decomposing goi = foiix) + ^i-^ii(^) in t^*^ lemma and absorbing 

F to G. If / is defined in an exterior domain we may assume go = by (i) above. 

Proof. We will prove the case n = 3. The proof for the general case is similar. By rescaling 
we may assume R ~ 1 and ||(a;)°/(.T)j|ioo = 1. We first consider the case supp / C -62- Choose 
a smooth ipit) : R ^ R, ip{t) = 1 for < > i?i = 3^^/^ and i;{t) = for t < The region 

{x : \xj \ < Ri} C Bi. Define 

f{xi,X2,X3)dx3, l3{xi,X2)= lim G3(x), (2.42) 

F3(x) = G3(.t) -7A(x3)/3(xi,a:2). (2.43) 
Then / = dsFs + ^' {x:i)h{x 1^X2) . Define 



G2(a;i,X2)= / h{xi,X2)dx2, hixi) = lim G2(a;i,a;2), (2.44) 

J-00 X2^+00 

F2{x) = ^'(X3)[G2(CC1,X2) - i^{x2)l2{xi)\. (2.45) 

Then / = dzF:i + ^2^2 + ^j' {X2W {x3)l2{x^). Define 



Gi(xi) = / ' /2(xi)da;i, /= lim Gi(xi) - / /, (2.46) 
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Fi{x) = i;'ix2)^'{x3)[Giixi) - iPix^)!]. (2.47) 

Then 

/ = ^3^3 + d2F2 + diFi + ,l:'{xiW{x2W{x^)I. (2.48) 
The last term vanishes if / / = 0. Also note Fj are at least as regular as /, 

\FAx)\ < II/IU-, and suppFj C {x : \xj\ < 2}. (2.49) 

For the general case, choose a smooth ip{x) supported in {a; £ K'^ : 2^^ < |.t| < 2} satisfying 
ip{x) > for 2-1 < |x| < 2 and E^-oo = 1 for x 0. Let (pkix) = (p{2-''x) for fc > and 

tP{x) = 1 - J2T=i ^{^-''x). Define 

afc = ^rj^//E^,, (fc>0). (2.50) 

Then I a/,, I < 2~^". Decompose 

oo 

f^^^fk, /o = /'0 + aoVo, fk = if + ak)ipk - ak-i^Pk-i, (fc > !)• (2.51) 

fe=0 

One verifies that 

' fo= I f, / /fe=0, (fc> 1), (2.52) 



supp/fe C {x : 2*^-2 < < 2'^^+!}, and \ fk{x)\ < Cixy". The rescaled function fk{x) = /fc(2'=a;) is 
supported in i?2. The previous case gives the existence of F^j^ j — 1,2,3, so that 

3 

h = Y.d,Fk„ \Fk,{x)\<C{x)-''+\ (2.53) 

Thus 

3 oo 

/ = /o+E^^-^:" P^-T.^^^' \F,{^)\< E 2-(-i)'^- < (2.54) 

J = l fc=l fc > ln|a:| 

This shows the lemma. □ 
2.5 Stokes system 

In this subsection we proved decay estimates for the Stokes system in E x K.'^ . Recall the fundamental 
solutions of the Stokes system (see [55] and [531 page 27]) 



S,(t,x) . r(,.).., + I ^^dy, Q,(t,.) = M^, (2.55) 



where T is the fundamental solution of the heat equation. It is known in |33l Theorem 1] that the 
tensor 5" = (Sij) satisfies the following estimates: 

\Did^S{t,x)\<CkA\^\ + Vt)-'-'-^'', (^,fc>0), (2.56) 

where indicates £-th order derivatives with respect to the variable x. 
A solution of the non-stationary Stokes system in M'^, 

dtW-Aw + Vp = g + VG, divty^O, {{t,x) (2.57) 

if g and G have sufficient decay, is given by 

w^A{g) + Q{G), (2.58) 

where 

{kg),{t,x)^l I Sij{s,x-y)gj{y,t- s)dyds. (2.59) 
Jo JK3 

{QGUt,x)^- [ [ dkS,,{s,x-y)Gjkiy,t-s)dyds. (2.60) 

Jo 
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Lemma 2.7 (i) Suppose G e L°°Xa+i, < a < 2. Then QG E L°°Xa with \\QG\\l--x^ < 



CallGII 



+ 1 ■ 



(a) Suppose g E L°°Xa+2, 1 < a. Then Kg G L°°Xi with \\kg\\L-=-Xt < Ca\\g\\L°°x^+2- 

(Hi) Suppose g G L°°Xa^2, ^ < ct, g is time periodic of period T > 0, and Jq J gdxdt = 0. Then 
Kg e L^X2 with \\Kg\\L^X2 < Ca||.9||L~x„+2 • 

Note that time periodicity is only assumed in (iii). 

Proof, (i) For QG defined by (|2.60p . using the estimate (|2.56p and integrating in time, we 
obtain 

\eG{t,x)\< / ^^{y)-^-^dyds < Ux), (2.61) 

Jo Jk3 (|a; -y\ + Vsr 

where 

1 



I : ^{y) dy, («>0). (2.62) 



-[V) 

\x - y\ 

If |a;| < 10, then /„(.t) < 1. If \x\ > 10, then 



(2.63) 



if < a < 2. We conclude /q(.t) < Ca{x) " and the estimate for QG. 

(ii) Due to Lemma 12.61 g can be decomposed as g = go + VG, where supp go G ^i(O) and 
||C?||loox„^i < ||5||l°°js:q+2 ■ Since VG can be treated as in the case (i), we consider the only the case 
that go ^ and G = 0. We may assume ||(7o||l~x„+2 — 1- By Young's convolution inequality 



\^g{t,x)\< I ||5y(s, •)|lL2+L~|l,gj(-,i- s)||L2nLic^s 
i{s~'/\s-^/^)ds<l. 



< 



(2.64) 



mm 







For \x\ > 2 we have 



(2.65) 



\^git,x)\< / / \S^j{x-y,s)\\gj{y,t- s)\dyds 

Jo J\y\<l 

poo n poo ^ ^ 

(iii) Continue part (ii) and assume \x\ > 2. Using J g{t,y)dydt = 0, we have 

{Kg),{t,x)= / {Sij{s,x - y) - Sij{s,x)) gj{y,t- s)dyds. (2.66) 

Jo JR3 

Using the mean- value formula and (j2.56p . 

\(Kg){t,x)\< f f \VyS^,{s,x-ey)-y\\g,{t-s,y)\dyds 

Jo J\y\<l 

<C / 71 ^dyds < G / ^ds < — 7, 

" Jo J\v\<ii\x-oy\ + v^r ~ Jo + ~ ixf 



(2.67) 



where 9 = 6{x,y,s) G [0,1], and we have used that gj is bounded and supported in \y\ < 1 and 
> 2. This shows ||A,g||Loojf2 < Ga\\g\\L°°x„+2 ^^'^ completes the proof □ 

3 Existence of flows with non-decaying boundary data 

In this section we prove Theorem ll.il By Lemma 12.61 we may assume fa = by absorbing /o into 
F. 
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3.1 Construction of a mild solution 

Proof of part (i). We first consider tlie case / = K. and denote 

A-* =i°°(M,i«'°°). (3.1) 

Let Ffc, k = 1, . . . , M ^ be the connected components of dVL, and clioose Xk in tlie bounded open 
region enclosed by Ffc. Let Ho{x) = j^. Note 

/ VHaix~xi)-N ^AirSki, Vfc,Z. (3.2) 

For each i G /, let 

H{t,x)=Y,9k{t)Hoix-Xk), gk{t) = ^ u^t) ■ N. (3.3) 
Note H is liarmonic in x e fl. Denote <(t) = u^{t) - \/H{t)\gn. By and 

<(t) • iV 0, Vfc, yt. (3.4) 



By Lemma [2.41 we can define Ei — £{u'l{t)) S C^(fi) with compact support in Br-^. 
Decompose 

It = « + £;, E = VH + Ei, p = Tr-dtH. (3.5) 
Then we have E\qq = u*, and w, tt satisfy 

- Aw + Vtt = 5(1;), divw^O, wlan = 0, (3.6) 

where 

giv) ■[F-{E + v)®{E + v)] - dtEi + AEi (3.7) 

and E satisfies ||i?||^3/2.co < ||u*||(72(gj2). Moreover, since dtEi is compactly supported in fi, we can 
apply Lemma |2. 61 so that dtEi can be written as 

dtEi = VFi with \\Fi\\l..^ < UrdtEiU^ < \\dtu4cHdn) 

for q > S/2 and a > 3. Therefore we can rewrite 

g{v) = VG(«) with ||G(«)||;,3/. < \\F\\xs,2 + (||«*!| + M^sf + \\u4. 

Here, = (72(3^^)), and we do not distinguish V and V- , since the difference does 

not play any roles in our argument. 

Following [37], we consider the fixed point problem 

v = ^v, {^v){t):=f e-(*-'')^PVG'(w)(s)ds, (3.8) 

— 00 

where the integral is defined weakly in the sense of Lemma 12.31 
Applying Lemma 12.31 we have 

ll'f^lU^ < \\G{v)\Us^. < IIFIU3/2 + {\\u4 + \\v\Usf + \\u4. (3.9) 



Similarly, 



Thus, if 



- $«|U3 < (||u,|| + \\v\\xs + Mx^)\\v - SlUa. (3.10) 



|F||;,3/2 + ||u,i| (3.11) 
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is sufSciently small, there is a unique fixed point of $ in the class 

Mx^<e, (3.12) 



and V G SC^(R, Ll-°°) by Lemma [ 

The case / = (0,oo) is proved similarly: We define gk{t) and Ei{t,x) as above for t > 0, and 
decompose u = v + E as m p.Sp . The vector field v{t, x) satisfies the initial condition 

v{0, x) = vo{x) uo{x) - E{0, x). (3.13) 

It is a fixed point of $ where 

= e-'-^^o + / e-(*-"'^PVG(i;)(s)ds. (3.14) 



Denote X"^ = L°°{M.-^, L'^'°°). By similar estimates, there is a unique fixed point of $ in the class 

\\v\\xi<Ce, £ = !|«o||l3.- + ||F|L.V2 + (3.15) 

if e is sufficiently small, and v G i3C^([0, od), i^^°°) by Lemma 
3.2 Equivalence to a very weak solution 

We now verify that our solution u ^ v + E is the unique very weak solution in the class u G 
BC^{I,L^'°°{n)) of (|TTT|) with the given data uq, u*, and /, with / R or / = (0,oo). 

A mild solution is a very weak solution. We first show that our solution is a very weak 
solution. Clearly u satisfies (|2.22p and divu = 0. It suffices to show (|2.2ip . 
By divergence theorem, for w G Vj, defined in (j2.20p . and fixed t, 

(u* , TV • Vw)f)^ ^{E, N- Vw)g^ = {E, Aw) + {VE, Vw) = {E, Aw) - {AE, w). (3.16) 

Thus, upon writing u ^ E + v, (|2.2ip is equivalent to 

(v , wt + Aw)dt = {vo , w{0)) - / (Fa , Vw)dt, Vw € Vj (3.17) 



for F2 = F2{v) = F - {E + v) <g) {E + v) + V Ei - Fi. RecaU VFi = dtEi. Above wo = if / = R, 
and vq is given by (|3.13D if / = (0, 00). 
Let 

V-{^G C2(n;R3),div7A = 0, Vlao = 0}. (3.18) 

Note VV'laa may not be zero for G V, and 2? is a dense subset of D{A). Choosing w{t, x) — 9{t)^{x) 
with 0{t) G C^{I) and ip eV, f^TU} implies 

iv , i;)e'dt = {vo , ^)6i(0) + J^e [{v, - (F2 , VV')] dt. (3.19) 

In turn this also implies (j2.2ip since linear combinations of such w{t,x) = 6{t)il;{x) with 9 G C^{I) 
and ■!/; G I? is dense in Vj in the norm J2a<i b<2 \\'^t'^\'^\\c^' ■ 

Denote = inf /. Being a mild solution v = oi (|3.8|) means ?j(i) satisfies 

(«(t),^) = (e-*^«o,V')+ A-F2(s),Ve-(*-^'^V)ds, G L^^^f^), Vt G /. (3.20) 



Plug this in the left side of For V e 2? and 6'(i) G C2((to, 00)), since [e'{t) - £-^{e{t) - 6l(t - 

e))] — > as e 0+ uniformly in t, 

-{v,il})e'dt= lim / {v{t) , il:)[-e{t) + e{t - e)]dt = (3.21) 

e-i-0+ 

= lim £-1 / (w(t + £) - v{t) , ij)9{t)dt = 4 + /4 (3.22) 
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where 



(3.23) 



+ y" 9{t) j\-F2{s),Ve-'^'-'^'^e-\e-'^ -l)ij)dsdt, 

He = e-^ J 6{t) y ^ (-^2(5), Ve-(*+^-^)^V)rfs (3-24) 

For t/- G X> C Z^CA), -A?/' = PAtA e L^^^'^ and 

||e-i(e-^^ - 1)7/. + AiZ-ILa/s,! ^ as e ^ 0+. (3.25) 
Thus, by LemmaOand v g BC„(/; Ll'"^), 

le^ f e{t){e-'\o,PAi^)dt+ I 9{t) f {-F2{s)ye-^*-''>^PAi;)dsdt 

J J Jto (3 26) 

= J e{t){v{t),PAi:)dt^ J e{t){v{t),A^)dt, 

as e — ?> 0+. Moreover, with t = {t + e — s) and by Fubini theorem. 

He = e-^ J 0{t) {-F2{t + £ - t), \'e-^'^i;)dT dt (3.27) 

= -e-^ f {( e{t)F2{t + e-T)dt,Ve-^^ilj)dT. (3.28) 



Since / e{t)F2{t + e - T)dt = / 0(^ - e + T)F2{t)dt and ^(^ - £ + r) - ^(t) ^ as £ ^ uniformly 
in t and t G [0, e], we get 

He +^^^ (J G{t)F2{t)dt, Ve-^'^ij) dr^O (3.29) 
as £ — > 0+. Since e~'^^ip is continuous in i^^^'^, we get 

lim //e = -( / e{t)F2{t)dt, VV'). (3.30) 

The above shows (j3.19p for i/j G V and 0(t) e C^((to,oo)). By approximation, (j3.19|) is also vahd 
ioTipeV and 0(t) e C^{[to,oo)). 

A very weak solution is a mild solution. We next show that if m is a small very weak solution 
in the class BCnj{I, L^'°°), then v = u — £' is a mild solution. This implies their equivalence, and 
also the uniqueness of small very weak solutions in the above class. 

Let u be small very weak solution in the class BCw{I, L^'°°), then v = u — E satisfies p.l9p . 
Taking 0{t) e C^{I), ([XTO)) is the weak form of 

^(«(t), V) = (v, A^) - (F2, W). (3.31) 

If we now take ijj = ^>{t) = e"'*^"*-'"^?^ with t < ti and rj £ V, we get -^ipit) = — A-0(t) and 
^(?;(t),e-(*i-*)^77) = -(F2, Ve-(*i-*)'^77) weakly, that is 

{v{t),e-^*'-'^^T])9'{t)dt = - f {F2it), Ve-(*i-*)^r/)0(i) dt (3.32) 



for any e{t) G Ci((inf Take to e [mf /,ti) and 9(1) = + I) - (t){^-^) where < £ < 1, 

(j){t) G C^K), (l>{t) = 1 for t < and (j){t) for t > 1. Send £ 0+. We have 6l(<) lto<t<t^ 
and; by continuity of (?j(i), e~(*i~*'"^77), 



(«(ti), 77) - («(to), e"(*^-*°)^r;) = - / (i^2(i), Ve-(*^"*)^7?) dt. (3.33) 
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we get 



If / = M, wc send to -oo. Since v S L°°(M; L^^°°) and e^^*!"*")-^?? ^ in L^^^'^ as to -oo, 

= - / ' (F2(t),Ve-(*i-*)^r;)dt, Vti G /, Vr; e P. (3.34) 

— oo 

If / = (0, oo), we take to — > 0+ and get 

{v{ti), ri) - {vo, e-'^^ri) = - / {F2{t),V e-'^'^-'^^ii) dt, Vti £ /, V?? e P. (3.35) 

Jo 

In either case v{t) is a mild solution and is the unique one we constructed in the previous subsection. 
3.3 Periodicity and spatial decay 

Proof of part (ii). Since u{t + T, a;) is another solution of (|l.ip - (|1.2p with the same data and 
estimates, we have u{t + T, x) = u{t, x) by the uniqueness of part (i). 

Proof of part (in). We now assume / = R and the stronger assumption of (iii). 
We first prove some a priori bounds. By Lemma 12.31 again for a fixed q € (3/2,3) and q*, we 
have 



Mx.' < \\Giv)\\x- 

< \\F\\x. + i\\E\\x^ + |K.|U3)(||iJ|U,. + \\v\\x,.) + \\u4. 

Hence, we obtain. 



(3.36) 

\\v\\x.' <\\F\U. + \\u4 (3.37) 



when e is sufficiently small. 

We now estimate w ^ vt. It satisfies 

wt- Aw + VTrt=VdtG, diviu^O, w\an = 0, (3.38) 

where 

dtG = {Ft - {Et +w)®{E + v)-{E + v)® [Et + w)) - dfPi + dtAEi. (3.39) 
For fixed 3/2 < q < 3 wc have 

WdtEWxo'nx^ + \\dfFi\\x,^x^/^ + W^t"^ Ei\\x,^j^3/2 < \\dtu4w^---iR,cHdn))- (3.40) 
Similar estimates show vt G n X'' and 

IktllA'SnAT'!* ^ £i ll-^llwi.^(R,L3/2nL<!) + ||u*||H'2.~(E,c2(aa))- (3-41) 

From now on we choose q = 2 and thus q* = 6 for convenience. 
Rewrite equation p.6p for w as a time-independent Stokes system, 

- Av + Vtt = g{v) - Vt, div w = 0. (3.42) 

We have v X^ f] Xi\ F <E X'^/"^ n A"', and -dtEi + AEi - vt £ X^ D X^' . By bootstrapping as 
in [16], also see [25], locally E and tt — (p{t) with a suitable (/)(t) are bounded uniformly in t. Let 

7^c«t {x e . < < + 1}. (3_43) 
We now replace 7r(i, x) by 7r(i, z) — (^(t) with a suitable (^(t) so that 

\vit,x)\ + \Vv{t,x)\ + \TT{t,x)\ <Cei, VteM, Vx e 7^™t. (3.44) 
Back to u = E + V and p = tt — dtH, we have 

|klU3n;t.* < Ce; (3.45) 

|u(<,a;)| + |Vu(t,a:)| + |p(t,x)| < Ce, Vt G M, Vx G 7^c„t. (3.46) 

The pointwise decay estimate of (iii) follows from the following lemma. □ 
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Lemma 3.1 (Spatial decay of time-entire solutions) For any Ri > 0, there are £i > and 
C > such that the following holds. Suppose u,p is a solution of p. II) with force f = fo + V ■ F 
for Ri < \x\ < oo and <: e R, and satisfies 

e :=sup ( sup {\x\^+"\fo{t,x)\ + \x\^+"\F{t,x)\} + \\u{t, •)||l3,oo(|,|>«^) + 

I N>«-^ (3.47) 

+ sup {\u{t,x)\ + \Vu{t,x)\ + \p{t,x)\}\ <si. 

i?i<|2;|<_Ri + l J 

We do not assume any boundary condition at \x\ = Then 

\u{t,x)\ <Ce\x\-\ {\x\> Ri; teR). (3.48) 

Proof. We perform a cut-off and the extend the solution for |a::| > i?i to entire M"^. Fix a smooth 
fimction C,{x) which is 1 for \x\ > Ri + 0.9 and for |a;| < i?i + 0.1. Let 

u = (u + \7r], divtt = 0, p~(p — dtri, (3.49) 

where, for each t, rj is defined by the Newtonian potential so that divu = 0, i.e., 

ri{t,x)=j -A,^..VC. (3.50) 
"'Ke„t 47r|a; - y| 

Note that u^p are defined for x E and satisfy the Stokes system for {t,x) e M^+^ 

dtu - Au + Vp = C/o + V • F3(ft) + /4, div u = 0, (3.51) 

wher(3 

Fsin) = - (u - Vr/) «) (u - Vr/), 

/4 = -(VC) • + V(^* • VC) + V • [(C' - On <E> u)] (3.52) 
- 2(VC • V)u - uAC + {u ■ VC)w + pVC. 

Note contains ii and global source terms, while f4{t,x) contains only local source terms, 

supp/4 C M X TZcuu IIMIl,- < Cei. (3.53) 

Consider now the fixed point problem w ~ for the map <f> from the class of vector fields defined 
on M X M'' into itself, defined by 

i^w) = A(C/o + fi) + e{F3{w)). (3.54) 

Here A and 6 are defined by p.59|) and (|2.60l) . respectively We want to show it is a contraction 
mapping in the class of small vector fields in L°°Xi. By Lemma [^771 

W^wWl^X^ < C\\Cfo + .fi\\L^X,+s+C\\F3\\L^x, 

< Cei + C(ei + \\w\\L^x^)\\w\\L^x^ (3.55) 

W^W - < C{ei + \\iv\\l=-Xi + \\w\\l'^Xi)\\w - 

Thus there is a constant Ci so that $ is a contraction mapping in the class 

{u;(i,a;) : R X R3 ^ K^, ||w||locXi < CieJ. (3.56) 

Thus there is a unique solution of w = in L°°Xi. 

By uniqueness of small solutions in L°°(R, L'^'°°(R'^)), we have w = u = C,u + Vr]. Thus 

|u(<,x)| < |w(t,x)| + |V7/(t,x)| < Ci£t^, (|a;| > i?i + 1). (3.57) 

□ 



^Eq. 113.5211 seems complicated because it is used for the construction of u without assuming decay estimates of u, 
and hence cannot allow u in _F3. 
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4 Spatial asymptotics of time periodic solutions 

In this section we prove Theorem 1 1.2 1 We start with a T-periodic solution {u,p) of (|l.ip - (|1.2p with 
force F for \x\ > R, satisfying the estimates assumed in Theorem 11.21 Proceed as in the proof of 
Lemma |3. II and let Ri = R, we fix a cut-off function ^(x) with VC supported in TZcut and define 

u ^ Cu + \7r], p = p( — dt'q, — A?7 = • u, (4.1) 

as in p.49p and p.50p . Then u and p satisfies p.5ip - p.53p . By Lemma [3.11 we have a;)| < 
Ce|a;|"i for |x| > R. 

Let b be the constant vector defined by (jl.l4p . 

b, = hm ^ / / T,j{u,p,F)n,dS^dt, (n, = i = 1,2,3), (4.2) 

p^oo T Jo \x\ 

where Tij{u,p, F) = pSij + UiUj — diUj — djUi — Fij. By the fast spatial decay of Vt?, the time 
periodicity, and divergence theorem, 



bj ^ lim / / T^j{u,p, F)n^dSx dt = ^ / / (C/o 

p^oo 1 Jo J\x\=p 1 Jo J\x\<Ri + l 



+ fi)jdxdt. (4.3) 



Also let be the corresponding Landau solution, and let 

U = C,U^ + U, P = CP\ (4.4) 
where U has compact support in TZcut, div U = —U ■ V^, and 

< C||C/.VC||M/-.a(7^,„,) < C\b\ < Cei. (4.5) 

Note that U, P satisfies 

- AC/ + V.((7®;7) + VP = /j^, div;7==0, (4.6) 

for a; € M'^, where 



f^=-AU~ 2WU'' ■ VC - U''AC + PVC + (C^ - C)V.(?7'' ® U'') 



Moreover, 



ifu)jdx^ I d,T,j{U,P)dx^ / T,j{U,P)^dSx 

'\x\<Ri+2 J\x\=Ri+2 Fl 



Let 



T 



JR'^ 



(4.7) 



+ V ■ (U ®U ~ CU" (E) cu") + {u" ■ vou" 

Wfifh^ < Cei, supp/^ c TZcut- (4.8) 



(4.9) 



/ T,,{U\p')^dSx = b,. 

J\x\=Ri+2 Fl 

W^U-C/, TT=p-P. (4.10) 



They satisfy, for (t,x) e 

atw- Au + V7r = .g + VG(w), divi; = 0, \v{t, x)\ < Ce2{xy^ (4.11) 
where (note G'(w) = ^3(14) + C7 t/) 

■'' = ^^° + ^^-^- , . (4.12) 

G(i;) = CP - (w - V?7) ® (w + C/ - V77) - [/ ® (w - Vr;). 

Because J^:i{fjj)idx = 6^, 

1 /-^ /• 

g{t,x)dxdt ^0. (4.13) 



18 



Consider now the fixed point problem v ~ $u for the map <f> from the class of vector fields defined 
on M X into itself, defined by 

$w = A.g + eG(i;). (4.14) 

Here A and 8 are defined by (j2.59p and (|2.60p . respectively. We want to show it is a contraction 
mapping in the class of small vector fields in L°°Xa, 1 < a < 2. By Lemma [2.71 



(4.15) 



||*^;||lo=x„ < C\\g\\L^x,^,+C\\Giv)\\L^x^^, 

< Csi + C||CF|U„ + C\\v - VyyllxJIk - VvWx^ + \\U\\x,) 

< Cei + C{ei + \\v\\l^xJ\\v\\l^x^ 
- ^v\\l'^x^ < C{si + WvWl'^Xo, + \\v\\l-^xJ\\v - v\\l'^x^- 

Thus there is a constant C2 so that $ is a contraction mapping in the class 

{v{t, x) -.RxR^ ^ M'\ \\v\\l^x^ < Cad}. (4.16) 

Thus there is a unique solution of w = $w in L°°Xa-, which agrees with it — J7 by uniqueness of small 
solutions in L°°{R, L^^°^{M.^)). 

Since v = (u + V77 — CU'^ — U, for |a;| > i?i + 1, we have 

\u{t,x)+VT]it,x)-U\x)\ < C2e{xy. (4.17) 

Since |V?7(t,x)| < C{x) ^, we have proven Theorem 11.21 □ 



5 Perturbed Navier-Stokes flows 

This section prepares a few lemmas for the proof of Theorem 11.31 in S|6l We first consider the 
solvability of the perturbed Navier-Stokes system (|1.19p in Proposition 15.11 and Lemma 15.21 We 
then prove a few estimates in Lemmas I5.3H5.4I 

We now recall a few notions related to self-similar solutions. The Navier-Stokes equations (|l.ip 
in / X J7 = (0,00) X with zero force enjoy the sealing property that if u{t,x) is a solution, then 
so are u\{t,x) := Xu{X^t, Xx) for any A > 0. li u\ = u for all A > 1, then the solution is called 
forward self-similar (SS), and u{t,x) = t^^/^it(l, t^^/^x). If ux = u for a particular A > 1, then 
the solution is called forward discretely self-similar (DSS), and it is completely decided by its values 
when T < t < TX^ for any T > 0. The existence of both type of solutions for small initial data uq 
with uq{x) = Aug (Ax) for all A > 1 or a particular A > 1, (also called SS or DSS), follows from Giga 
and Miyakawa [T3], see also Cannone, Meyer and Planchon [5] and Cannone and Planchon [B]. 

We now state our result on the solvability of the perturbed Navier-Stokes system which is useful 
to describe the time-asymptotics of solutions for the initial-boundary value problem (jl.l[) - (jl.3p in 
the next section. In particular, it implies the existence of the self-similar solution in Theorem 11.31 



Proposition 5.1 (Perturbed system) For any < 77 < 1 there is Eq = eo{r]) > such that the 
following holds. Let U,U,Wo be vector fields in M.^ satisfying sup^ |x|(|C/(x)| -|- |f7(a;)| + |u'o(a;)|) < 
e < eo and divwQ ~ 0. Then there is a unique solution w{t, x) of the perturbed Navier-Stoke system 

dtw ~ Aw + V{w(^ w + U (^w -\-w'S>U) -\-Vp ^ 0, divw = (5.1) 

in MP with initial data w{0) = Wq, satisfying 

\wit,x)\ < ei\x\ + Vty^+'^lxl-". (5.2) 

The pressure satisfies |ip(i)||L»,°=(R3) < e^t^^+A for any | < s < 
If furthermore sup^ \x\'^{\\7U{x)\ 4- |Vf/(a;)| + |V?i;o(a;)|) < e, then 

\Vw{t,x)\ < e{\x\+Vi)-^+'"\x\-^-'K (5.3) 

Comments for Proposition \5.1i 
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1. U and U need not be divergence- free or self-similar. If we assume in the following that both 
U and U arc SS, the solution set of (|5.1|) has the same scaling symmetry as the usual Navicr- 
Stokes flows: If w{t,x) is a solution of (|5.ip . then so is 

wxit^x) = Xw{Xt, A^x) (5.4) 

for any A > 0. By this scaling symmetry and the uniqueness part of Proposition 15.11 w is 
forward SS (or forward DSS) if wq is SS (or DSS). 

2. If w is forward self-similar, then w{t,x) = W{xl\ft)l\Jt with W{x) = ^(1,2;). In the case 
L/ = f/ = and is SS and small, we have W G L^^°° and expect |W(j/)| < , i-e., 

|w(t,a;)| < e(|x| + (5.5) 

3. For general initial data |wo(a^)| < e|a;|~^ and ?7 = f/ = 0, we can construct solution of (|5.ip 
satisfying the same bound (j5.5p . However, this is impossible if [/ or [/ is nonzero. See Remarks 
(i) and (ii) after the proof of Proposition 15.11 

4. In the case U and U are nonzero, and all C/, U and are SS and small, Wiy) satisfies 
div T4^ = and 

\^ ■ - + '^(^ ®W ®W + W ®U) + ^P (5.6) 

Since W £ L^'°°, we expect |W(2/)| ^ for large y. Due to the local singularity of U and 

U, and the local analysis of [IS], we can show |M^(y)| < ^1^1'' for small y and some 77 > 0, 
with a smaller e needed if we want a smaller 77. Thus 

\W{y)\<e{\yr + \y\)-\ (5.7) 

which corresponds to (j5.2p for w. For general data wq G Xi, Proposition l5 . 1 1 asserts the unique 
existence of solutions w satisfying the bound (|5.2I) . 

In order to prove Proposition 15.11 we start with a lemma. 

Lemma 5.2 Let n e N, b > 0, c > 0, b + c < n, ^j. > 0, X > 0, and t>0. Then 

I {\x-y\ + X)-'•\x-y\~'{\y\ + Vty■"-^'dy^V^\\x\+X + Vt)-\x + Vt)-''. (5.8) 

Proof. We may assume t = 1 since the general case follows from the change of variables 
y = ^/iy, X = \fix and A = \/tX. Denote the integral with < = 1 as Jix) and p = |a;| -I- A. We want 
to show J{x) ^ {p)^^ {x)^^ . Let a = /i + 6 + c. If p < 10, we have 

Jix)< f \x-y\-'>--dy+ [ \y\-"-^dy<l. (5.9) 

J\y\<20 "'|y|>20 

We also have J{x) > /|y|>2o l?^! " "^2/ ^ 1- Assume now p > 10. We have 

J{x) = [ + [ =:J, + J2. (5.10) 

J\y\<p/i J\y\>p/'i 

For J2, by replacing the factor {\y\ + 1) of the integrand by |j/| and rescaling y — > y/p, 

J2<p-"/' {Ix^yl + Xr^lx-ym-^-^+'+^dy, (5.11) 
J\y\>i/i 

where x = xj p^ X ~ X/ p and |a;| + A = 1. The integral is of order 1 by the previous case p < 10. 
Thus J2 < 

For Ji, because \x ~ y\ + X> p/4 when \y\ < p/4, 

■h<p-'l \x-yr{y)-''-^''-'^^'dy. (5.12) 

JR" 
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Repeating the previous argument with A = and (6, a) replaced by (0, a— b), we get Ji < p ^{x) . 
We also have Ji > P"''/|y|<5/2 \^^v\~'''^y ^ P^^{x)~''- 

Since p-" < p-^xy", we get J{x) - p^^xy". □ 

Proof of Proposition [STTl 

Equation (jS.ip can be written in the integral form using the heat kernel F and the Stokes tensor 
(|235ll . 

w = wl + ujn{w), (5.13) 



where 



WL = e Wo, WL{t,x)^ T(t,y)wo^i{x-y)dy, 



WNiw){t,x) 







dkSij{s, y)Fkj{x ~y,t- s)dyds, 



(5.14) 
(5.15) 



and = F(w;) = wi^w + Ui^w + wi^U. Since e ^' < Ck (p) ' for any > 0, we can take k — A 
and have 

r(<,y) < i"'/'(x/%/<)" =Vti\x\ + Vi)-^. (5.16) 

Suppose |wo(a;)| < e(|a;| + A)^\ A = 0, 1. We have 

\wLit,x)\< j Vt{\y\ + Vt)-^e{\x-y\ + \)-Hy. (5.17) 

By Lemma [5. 21 with (n, /i, &, c) — (3, 1, 1, 0), 

\wL{t,x)\<e{\x\+X + Vt)-^. (5.18) 
Suppose further |Vti;o(a;)| < e(|x| + A)^^. By Lemma [Ol again with (n, 6, c) = (3, 1, 2, 0), 

|Vw;L(t,x)|< y"^/^(|2/| + ^/^)-4e(|:E-y| + A)-2^^J/ <e(|a;|+A + \/t)-2. (5.19) 
For < < 1, define two norms for functions on 



j3. 



\I\\yi = sup 



{\x\+VtY-'^\x\^\I{t,x)\ 



ll/b. = ll/b.+ sup \{\x\ + ^f-^^\x\'+^\Vf{t,x)\ 



(5.20) 



(5.21) 



Estimates (jS.lSp and (|5.19D show < Cie (resp. H^Lba — Cie) if ||wo||xi < e (resp. |lwo||xr 

||Vwo||x2 < s) for some Ci. 

We now estimate the nonlinear term wn{w). We will show 



\\wn{w) - WN{w)\\yi < e\\w - w\\y-^ 
if II w II < e and ||w)||j;i < e, and 



\'>^N{w)\\y., < e\\w\\y^ 



(5.22) 



(5.23) 



if 



M\y2 



< e. Note (|5.22p implies ||u'Ar(w)||;yi < ellu'llj;^ by taking w — Q. These two estimates 



imply that the map 



w —J- Wl + wn{w) 



(5.24) 



is a contraction mapping in the class of vector fields defined on R+ x K'^ with H^bi < 2Cie 
(resp. ||u'||3;2 < 2Cie) if ||wo||a'i < £ (resp. ||wo||xi + ||Vwo||a'2 < £) and e is sufficiently small. 

If \\w\\y-^ < £, ||w||3;i < £, and ||C/||xi + ll'^llxi < £, then 



|F(u;) - F{w)\{t,x) <e\\w- w\\y,{\x\ + Vt)-^+^ \x\-^-^ . 
If \\w\\y, < e and \\U\\x, + \\U\\x, + \\yU\\x, + \\yU\\x, < £, then 

\VF{w)\it,x) < e^i\x\ + Vt)-^+'"\x\-^'\ 



(5.25) 



(5.26) 
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Thus, using the definition of wn and V5-cstimate (j2.56p . to prove (|5.22p and (|5.23[) it suffices to 
show 







I + ^/^)-Wx -y\ + Vt^r'+'' \x - yr'^"" dyds < {\x\ + Vi)-^+n Ixl'"^'"" (5.27) 



\y 

for fc = 1, 2. Decompose the integral as 



ft 

I + 11 :^ I I ... dyds + I I ■■■ dyds. (5.28) 




We first estimate /. Since (jx — y| + \/t — s) > t, 

I < Vf'^'' r' I {\y\ + V^)-^ \x ~ y\-'-'' dyds. (5.29) 
Jo JK3 

By Lemma 15.21 with (ri, /i, 6, c) = (3, 1, 0, fc + i]) and t replaced by s, 

I<V-f'^'['^s-'{\x\ + V^)-'-'ds^J^['^—-^^^—. (5.30) 

where we used the scaling s = t. Using y7(i_|^^)fc+^ ^ for fc + 7y > 1, (thus we need 

77 > for fc = 1), we get 

1 

/ < , ^ (5.31) 

N'+^-^N + Vt) 

which is bounded by (|a;| + \/t)^^^^' \x\^~'^~'' . 

Next we estimate //. Bounding the factor {\y\ + \/s)~^ by {\y\ + Vt)~'^ and then integrating in 
time, 



II < / {\y\ + V~t)-%\x-y\ + Vmy^"-\^-y\Y^']\^-y\''^'dy. (5.32) 
Using {A + B)^+''' - ^1+'' - (A + S)''B - A^'B + 5^+'' for A, B, > 0, 

/ (|jy| + \/t)^^[k-y|)-^-Vt+|x-yr'^-"Vt'^"]d2/. (5.33) 



By Lemma 15.21 with (n, a, &) = (3, 1, 0) and c = fc or c = fc + ry, 

^^<(Nl+v<)-^ (5.34) 

which is also bounded by (|a:| + \/t)^^^^ \x\^^'^^^ . Summing up, we have shown (|5.27p and thus the 
existence of w satisfying (|5.2p and (|5.3p . 

It remains to show the estimate of the pressure, which follows from its equation 

~ Ap djd^iwiWj + U.,wj + WiUj) in R^ (5.35) 

the Calderon-Zygmund estimates, and F(w)-estimate in (|5.25p with w = 0. □ 

Remarks, (i) If t/ = [/ = and |wo(-'J;)| < £{\x\ + A)^^ with A = or A = 1, then one can 
construct solutions in the class 

\w{t,x)\ < e{\x\+ X + Vi)~\ (5.36) 

Indeed, in this case, we have | J^(t, x)| ^ e^(|2:| + A + Vt)^'^ and can bound WN{'w){t, x) by Ce^(|x| + 
A + \/t)~^ using 



Ji\y\ + V^yWx - 2/1 + A + Vt~s)-^dyds < Vt{\x\ +X + Vt)- 



(5.37) 
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(ii) If J7 or J7 is nonzero, we need to take 77 > in (j5.2p . If r/ = 0, we have |i^(t,a;)| < + 
and need 

f f i\y\+V^)-Wx~y\ + Vt^)-'\x-y\~'dyds<{x + Vt)-\ (5.38) 

Jo JR3 

However, when |a;| <C Vt, 

I> f'^ I {\y\ + V~srWt)-^\x~y\-^dyds. (5.39) 

Integrating in ds first, 

/> / \y\-^{Vt)-'\x-y\-Hy. (5.40) 

Restricted in the subregion 2\x\ < \y\ < y/t, 

I> I \yrHVt)-'dy^[Vt)-'\og^ (5.41) 

It is larger than {x + \/t)^^ by a factor log when |a;| ^ \/t. 

(iii) Since V^F ^ L\, and V^S* ^ LiJ, ^ j^^, there is no suitable integral formula for V^w, and 
the above method does not allow us to estimate \/^w pointwisc. 

(iv) If all U, U and wq are self-similar, then w{t,x) = W{x/^/t)/^/t with W{y) satisfying the 
elliptic equation (|5.6p . and one can estimate higher derivatives of W. It is not clear if wq is DSS. 

In the rest of this section, we give two lemmas to be used in the next section. 

Lemma 5.3 Let TZcut o,nd cut-off function Q be as in section\^ There is a linear map 

A: {we Ll^{R^;R^), divw = 0} ^ {w e L,io^(R^ M^), suppw C TZcut} (5.42) 

such that, for any 1 < g < 00 and q <r < 00, w = Aw satisfies divw = —VC • w and 

l|w||H'i>5-'-(K,„t) + I|VPw||l,,.(o) + ||Pw||L3/2,oo(o)nL5i.-~(0) < C'9l|w||L5.-(K„,t) (5-43) 

where P is the Helmholtz projection on L'^'^{Q,;R^), ~ q* if q < 3, and qi, ~ 100 if q>'S. 

Proof. The usual construction of the solution of the problem divu = / (see [TOl Thm. III. 3.1]) 
gives such a linear map with ||it'||vKi.9(7?,c„t) < ^qll'^WLUKaut)- ^^^'^ bounds in Lorentz spaces follow 
from interpolation. It remains to show ||Pw||iyi>9''-(o) £ C!q\\w\\Lq.T(^-ji^^_^y Decompose w = Pw + Vtt 
where the scalar function tt solves the following Neumann boundary value elliptic problem: 

Att = -VC -w in n, Vtt • TV = on dQ. (5.44) 

Due to LP theory for elliptic equations with Neumann boundary condition, (by partition of unity 
and boundary estimates in [1], also see [32]), we have 

l|V^7r||i,,r(j2) < ||VC • < l|w||L9.'-CRe„t)- (5-45) 

Combining estimates, 

||VPw|U,,,.(o) < \\w\\w^.-,.rin) + < \\w\\L,.r(n^^,). (5.46) 

To bound |lP?i'||L3/2,oc(Q)Q2:,g6,oo(Qj amounts to bounding Vtt is in the same space. This can be 

shown by cut-off: u = ttC is in W^^^'^ , Cf^^ outside of TZcut, and solves Am = / for some / G L''-^ with 
compact support. Thus |Vm(x)| < C|a;|~'^ for |a;| large using Newtonian potential. This completes 
the proof. □ 
Note: Although / = Am is a divergence, we do not know if J f = (which would imply extra 
decay for m) since Am is not integrable and we cannot use divergence theorem. 
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Lemma 5.4 Let w be the solution in Provosition \5.1\ with wq and Vti'o satisfying the stated esti- 
mates. Let w ~ Aw be defined by Lemma \5. 3\ Then, for all 3 < q < qi, 



\\w{t)\\w^,^.^ + \\ e-^'-'^^Pdsw{s)ds\\L^..^ <Ce{ty (5.47) 



where = - ^ > 0. 

Proof. Denote fi = i-^. Note cr < ji ii q < qi. By Lemma and Proposition lS.ll 

ll^Wlk^-- < lkWllL--(7^.„.) < II^^WIIl^(7^.„.) < e(t>'^ (5.48) 
For the integral, using integration by parts, we have 

e-it-<>)Apd,w{s)ds = Pw{t) - e-'^^PwiO) - IC{t) (5.49) 

Jo 

where ^ 

/C(t) = / Ah-^'-'^'^AiPwis)ds. (5.50) 



Choose p so that fi — |(^ — i), with ^ = | ^ (We need q > 2{v-fi) ^ T+i;') Lemma 
and Proposition 15. 11 

WPwimL.-^ < WVPwim , ^.^ < e{t)-^. (5.51) 

||e-*^Pu;(0)|U„.^ < t-n|Pu;(0)||^^,^ < et-^, (5.52) 

while for i < 1, ||e-*^Pi&(0)||L.!-~ < ||Pw(0)||i,<! < e. 

To estimate /C(t), let Z = {ip e L^M ; < 1}, with l/q + 1/q' = 1 and r = (g')* i.e. 

1/r = l/q' - 1/3. We have 

i'^||/C(t)||L9.~ = sup(t^/C(t),(y5) ==sup(t^ /" (s-^A^e-(*-')'^¥5,s^A^Pi&(s))ds (5.53) 

ip Jo 



< 



supi^ /" s'^^'\\Ah~^'-'^'^ip\\Lr.lds sups^'\\A^Pw{s)\\^r'.^ (5.54) 

V Jo s<i 

<esupt'' f s-^P^e-(*-^'^(^||i,.,ids <e|l(^||^,M, (5.55) 



where we used (|2.13p . This completes the proof. □ 



6 Time asymptotics 

In this section we consider large time asymptotics of the solution close to the periodic solution. We 
obtain Theorem ll.3l as a consequence of a more general result about the asymptotics: 

Theorem 6.1 (Time asymptotics) For any T>Q, 6>0, rj>0 and 3 < qi < there is 

£2 > such that the following holds. Let and f be time-periodic data satisfying (|1.9p . Let uq £ Xi 
be initial data satisfying \\uo\\i^3.<x, < S2. Assume that there exists a vector field uq G Xi such that 

||"olUi<e2, Vwo e ^2 and e -.^ \\uo ~ uo\\ ^ 3 _^<e2, (6.1) 

then the solution u in Theorem ] 1.1\ (i) can be decomposed as the follows: 

u = Q + w + r. 

Here Q is the periodic solution for the data and f in Theorem \1.1\ (Hi). The term w is the 
unique solution of the perturbed Navier-Stokes system (|5.ip in M'^ with U = tj =U^ and initial data 
Wq = Uq — U^ in Theorem \5.1[ where is the Landau solution corresponding to Q given in Theorem 
with a = I + 6. The term r satisfies the following decay estimate: 
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Comments for Theorem \6.1[ 



1. Theorem 11.31 is a special case of this theorem when uq is self-similar (or (-l)-liomogeneous). 
We are also able to show that w is DSS when uq is DSS as we mentioned in the comments of 
Proposition 15.11 

2. The restriction g > 3 is related to the fact on the lack of the coercive estimate for the Stokes 
operator in LP'°°{^1) for large p in the exterior domain. See Lemma [2.11 Note that it is not 
the case when £7 = R'^ or when ft is bounded. 

3. By taking (5 + 77 small, we can choose qi < 00 arbitrarily large. But our method does not allow 
qi = 00. 

4. Our method allows us to consider u^, = + where is time periodic and decays in 
time. We assume = for simplicity. 

Proof. Suppose dil belongs to the region < i?o < |a;| < Ri and let TZcut and cut-off function 
C(.t) be defined as in section [31 To prove Theorem 16.11 we will study z = u — Q — w instead of 
r = u — Q — w where w is a cut-off of w which vanishes near dfl: 

v = u — Q = 'w + z, w = (^w + w, (6-2) 



where w with compact support in TZcut and div w — —w ■ \/( is given by Lemma 15.31 Since the 
difference w — w is localized near the boundary, it is easy to show that it decays in time. Compared 
to (j3.49p . this choice ensures z\q^ = but introduces dtw as a source term in dtz equation. 
Specifically, v and w satisfy 

9tw - Aw -t- Vpi = -VK,, divw = 0, (6.3) 

dtW - Aw + Vpo = -VFuj, divw^O, (6.4) 

with 

Fy = {v + Q)®v + v(»Q, = (w -f- C/'') w -t- u) (g) ^7^ (6.5) 

and 

v{0) = vo uo - g(0), vlan = 0, w{0) = wo := uo - U\ (6.6) 
Thus, w defined by (j6.2p satisfies the Stokes system for {t, x) E K+ x R'^ 

dtw-Aw + \/{(po)^-\7{CFu,+\7w)+f + dtw, divu; = 0, (6.7) 

with u)(0) ~ wq :~ ^wq + wq, where 

/ = (VC) • F„, - 2(VC • \7)w - {AOw + poVC. (6.8) 

The vector field z = v — w = u — Q — w, defined by (j6.2p . satisfies, for p^ = pi — C,pQ, 

dtz - Az + Vpa = VF^ - / - dtw, div z = 0, {x£n) (6.9) 

z(0) = zo, z\aa = Q, (6.10) 

where Fz — CFw — Fy + Ww and 

zo P{vo - ^0) = P[{uo - Ho) + {U'' - Q(0)) + ((1 - Qwq ~ wo)\ (6.11) 
and, using (|1.15l) with a = 1 + 5, 

llzoll 3 „ < e. (6.12) 

Note we can add the Helmholtz projection P in the definition of zq since z(0) = zg is understood in 
weak sense. Since 

Fy - Fy, ^z ® z + [w + Q) ® z + z ® [w + Q) + [Q - U^) ® w + w ® {Q - U^) 

, , (6.13) 
- {w-w)<^{w + U") ~{w + U") <»{w- w), 
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we can decompose as 

F, = z(E}z + Fi{z) + F2 + F3 + Fi (6.14) 

where 

Fi{z) = z <^ {w + Q) + {w + Q) <E) z, 
F2^w(E){Q- U^) + {Q- U'') (E) w, 
F3 ^ {w - w) {w + U^) - {w + U^) ®{w- w), 
F4 = {C- l)F^ + Ww. 



(6.15) 



Note that F3 and F4 have compact supports. 
We wih prove 

snpt'^\\z{t)\\Li.oc. < e for [8,(71], (6.16) 
t>o 

with 

a = a(g,,5) = ^(i-i) + ^>0. (6.17) 

By unique existence of m G L°°L^'°° fThcorcm ll.ip and Holder inequahty, it suffices to prove the 
existence of z satisfying (j6.9l) , (|6.10p and 

\\z\\z sup{t^/^z{t)\\L3.^ +e'\\z{t)\\L.i:^} < e (6.18) 
t>o 

where cti = cr(gi, S). 
Denote 

*(/)W= f e-^'-^^^Pf{s)ds (6.19) 
Jo 

where the integral is in the weak sense, 

i^Um^ip)^ f\f{s),e-^'-^'>\)ds, ycpeC^in). (6.20) 
Jo 

A solution z{t) of (j6.9p and (|6.10p is a fixed point of the nonlinear map $ defined for z G Z, 

$z(t) = z\t) + «'(V[z (g>z + Fi{z)]){t) (6.21) 

where 

z\t) = e-'^zo + ^(V[F2 + + Fi] - f-dtw){t). (6.22) 
We will show that for e sufficiently small, 

||$z||z<C£, \\^z-^2\\z <Ce\\z-S\\z, (6.23) 

which guarantees the unique existence of a fixed point z of $ with \\z\\z < Ce. 
By Lemma [221 (i): have 

< lko||^^.„. (6.24) 



Consider the terms of divergence form. By duality, we have 

I1*(V^^)|1l,- = sup \{^iVF),ip)\. 

By (jTTOl) with l/r = 1/q' - 1/3 and 1/r' = l/q+1/3, we have 

|(vI/(VF) , ^)| < /* \\F{s)\\^.,.^ II Ve-(*-^)^^|U.ids 
Jo 



(6.25) 



< fsupr||F(t)||^,.,„^) / s-n|Ve-(*-^)^^|U..id5 (6.26) 
\t>o J Jq 



< 



supr||F(t)||^.,,„ t-^l^ll^,,, 
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In the above we need 1 < g' < 3/2, thus 3 < q < oo. 
We now estimate F ^ F^. First, 

f'Wz ® z + < r |lz|l,,oo(|lz|l3,oc^ + ll^'lU.cc^ + llQ||3,oo) < e'. (6.27) 

Next, since w — tv = (1 — + w is compactly supported, it foUows from Lemma l5.3l and Proposition 
15.11 that, for 1 < g < oo, 

\\W - w\\L,.^i^n} < Wwil - C)llg,oo + \\lb\\q,oo 

£ - C)llg,oo + l|iy||L<!.~(TCe„0 (6.28) 

^ l|w'IU-(f2) < £{t)^^- 

Also note 

lh(i)llL'--(K3) <et-^+*, (3<r<3/r;). (6.29) 
With l/q\, = l/q — 6/3, (hence ^ — ^ = <^{q,S)), we have 



(6.30) 



Recall we take a = 1 + 5 in Theorem 11.21 so that Q — G L^+'''°°. We have also used that 
< a < ^-^^ for 3 < q < qi, which also implies 

llu; - w||l<j,oo(j2) < et"'^. (6.31) 

Thus 

i"||i^3||r'.co <t''\\w~ ^i)|U,oo(||^«||3,co + ||^i||3,cc^ + \\U%,oo) < ■ (6.32) 

Finally F4 — (C — 1)F„, + Ww can be estimated as follows: 

||-F4||r',oo < {\\w\\3,cc + 1 1 C^'' i 1 3,oo ) 1 1 (C - l)w'!lg,oo + {n^^t) 



(6.33) 



Next we consider ^'(/). Denote ri = 2+t+ri ^'^'^ ^2 = Note 1 < j'l < 7'2. By decay 

estimates, 

ll*(/)(t)IU,oo < (\t-s)--^\\~f{s)U^^ds, a2 = |(--i). (6.34) 

Note < (72 < 1 since q < qi < By Proposition 15. 1[ on 7?.cut we have |w| + \\7w\ < 
l^^l < +e|u;| < £2(t)-H^ and ||po(t)||r.,co < Thus, using supp/>) C TZcut, 

\\m\\r,,oo < \\fis)\Uoo < e{s)-^+^. (6.35) 

We get 

ll*(/)(t)IU.oo < f {t- sr--e{s)---+id.s < et---+i+i = et-r (6.36) 
Jo 

Finally the estimate '^{dtw) is proved by Lemma [5.41 

The above shows the first assertion ||$2||z < Ce in (j6.23p . The second estimate in (|6.23p on 
W^z — $z||z is proved similarly. This finishes the proof of Theorem 16. II □ 
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